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^ ' Abstract 

CD ■ We consider boson star solutions in a D- dimensional, asymptotically anti-de Sitter spacetime and 

investigate the influence of the cosmological term on their properties. We find that for D > 4 the boson 

■ star properties are close to those in four dimensions with a vanishing cosmological constant. A different 
£N) ' behavior is noticed for the solutions in the three dimensional case. We establish also the non-existence 

of static, spherically symmetric black holes with a harmonically time-dependent complex scalar field in 
any dimension greater than two. 

1 INTRODUCTION 

. Recently a tremendous amount of interest has been focused on anti-de Sitter (AdS) spacetime. This interest 

is mainly motivated by the proposed correspondence between physical effects associated with gravitating 
fields propagating in AdS spacetime and those of a conformal field theory on the boundary of AdS spacetime 

"o mm 

Being a maximally symmetric spacetime, AdS spacetime is also an excellent model to investigate questions 
of principle related to the quantisation of fields propagating on a curved background, the interaction with 
bi), the gravitational field and issues related to the lack of global hyperbolicity. 

In view of these developments, an examination of the classical solutions of gravitating fields in asymp- 
totically AdS (AAdS) space seems appropriate. 

Recently, some authors have discussed the properties of gravitating SU(2) nonabelian fields for A < 
[HQ] 03- They obtained some surprising results, which are strikingly different from the results found in the 
asymptotically flat case (for example the existence of regular dyon solutions without a Higgs field). Insofar 
as the no-hair theorem is concerned, it has been shown that there exist stable black hole solutions in SU(2) 
Einstein- Yang-Mills (EYM) theory that are AAdS OH]. However, the same system in the presence of a Higgs 
scalar presents solutions with very similar properties to the asymptotically- flat space counterparts H3 |S] • 
More recently it has been shown that the U(l) Higgs field equations have a vortex solution in four-dimensional 
AAdS 9 . Static spherically symmetric solutions of general relativity with A < 0, coupled to a real scalar 
field, have been studied in Surprisingly enough, the solutions have finite energy although there are 

naked singularities and present interesting properties in the context of the AdS/CFT correspondence. 

All of these studies assume that the matter fields have the same symmetries as the underlying spacetime; 
in particular, the scalar fields are assumed to be invariant under the action of the stationary Killing vector. 
However, the consideration of static configurations in which a complex scalar field varies harmonically with 
time, with a static energy-momentum tensor, yields new kinds of regular configurations, the so-called 'boson 
star' solutions. 

The study of boson stars can be traced back to the work of Kaup |TT| and Ruffini and Bonazzalo |12| 
more than 30 years ago. They found asymptotically flat, spherically symmetric equilibrium solutions of the 
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Einstein-Klein-Gordon equations. These boson stars are macroscopic quantum states and are only prevented 
from collapsing gravitationally by the Heisenberg uncertainty principle. Although boson stars have many 
similarities to their fermionic counterparts, there are many interesting differences. For example, boson stars 
also exhibit a critical mass and critical particle number. Later works considered the more complicated 
possibilities of the presence of a self interaction for the scalar field ^3 ^] E| or a nonminimal coupling of 
the scalar field to gravity [T^]. Jetzer and van der Bij extended the model to include a U(l) gauge charge 17 
(see also Boson stars in the presence of a dilaton or an axidilaton have also been studied by various 

authors [T^], as well as boson- fermion stars |2U| . 

All these models have demonstrated the same characteristic: new interactions tends to increase the 
critical values of mass and particle number, although the particular values are very model dependent. The 
stability against perturbations around the equilibrium state has been discussed also by a number of authors 
[2]- [Sj. An extensive review of the boson star properties is given in [22], |2fi| . 

All of these solutions assumed a spherically symmetric ansatz; a more recent development is the study 
of axially symmetric boson star configurations |27l 1281 Asymptotically flat, rotating boson stars have 
been obtained in |301 131j . 

The recent interest in these solitonic objects is largely due to the suggestion that the dark matter could 
be bosonic in nature. With scalar fields often used to model early universe, there is the possibility that such 
solutions could condense to form boson stars. 

Most of the works on boson star properties have been carried out on the assumption that spacetime is 
asymptotically flat. Because of the physical importance of these objects, it is worthwhile to study generaliza- 
tions in a different cosmological background. We may hope that similar to a nonabelian theory, the different 
asymptotic structure of the spacetime will affect the properties of the solutions leading to some new effects. 

The goal of this paper is to reexamine the basic properties of soliton stars in the presence of a negative 
cosmological constant. The solutions we are looking for are the asymptotically AdS analogues of the Ruffini 
and Bonazzalo configurations |12j . To our knowledge, a discussion of this case is not available in the literature. 

We find that the basic properties of a four dimensional soliton star will remain basically the same 
even in the presence of a negative cosmological constant. However a nonzero A term in the action implies a 
complicated power decay of the fields at infinity, rather than exponentially. Also, as expected, the parameter 
of the solutions found in Ref.|12| (mass, particle number, effective radius) remains no longer valid and new 
values are found for every value of A. The main effect of a nonzero cosmological constant is to decrease the 
maximal mass of the star. 

Motivated by recent interest in higher dimensional gravity with negative cosmological constant, we look 
also for spherically symmetric boson stars with A < in spacetimes with a number of dimensions D ^ 4. It 
is always of interest to see how the dimensionality of spacetime affects the physical consequences of a given 
theory. 

The higher-dimensional boson star solutions might be relevant in the context of superstring theory as 
well as for understanding how the behaviour depends on the dimensionality of the spacetime. In particular 
we are interested in the five-dimensional case. To our knowledge very little is known about the properties of 
these type of configurations in D =/= 4, although the exact boson star solutions in three-dimensional gravity 
with A < were found in |32| . A qualitative discussion of the boson star properties in the context of brane 
world models with compact large extra dimensions has been done in |33| . As argued there, the properties 
of large soliton stars will remain similar even in the context of the brane world and large extra dimensions 
scenarios. However, at distances smaller than the compactification radius, the properties of a boson stars 
are quite different than in the (3 + l)-dimensional case. 

We find that for D > 3 and a negative cosmological constant, the properties of the solutions are similar 
to the well-known four-dimensional asymptotically flat case. The total mass parameter M rises quickly with 
the initial value of the scalar field at the origin ipa to a maximum, then drops and has several local minima 
and maxima until it approaches an asymptotic value for large tpo- A similar behaviour is found for the 
particle number N. The solutions for a vanishing cosmological constant are discussed as a particular case, 
since they present the same general pattern. A separate section is dedicated to the rather special case D = 3. 
Exact solutions in this case (in the limit of large self interaction) are know to exist. However, the properties 
of the boson star solutions in AdS^ spacetime are different as compared to the higher dimensional case. We 
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notice the existence of a maximal allowed value for the central density and the absence of local extrema for 
N and M. Also, no solutions are found in the asymptotically flat limit. 

One may ask about the possible relevance of these macroscopic quantum states within the AdS/CFT 
correspondence. Since the proposal of Maldacena's conjecture, which gives a correspondence between a 
theory of gravity in AdSp spacetime and a field theory on its (D — l)-dimcnsional boundary, much intense 
work has been devoted to get a deeper understanding of its implications. The AdS/CFT correspondence for 
a scalar field has been discussed by a number of authors, however without considering these type of time- 
dependent solutions. We present in this paper only a preliminary approach to this question. The boundary 
stress tensor and the associated conserved charges are computed for three, four and five dimensions. In three 
and five dimensions, the counterterm prescription |34| gives an additional vacuum (Casimir) energy, which 
agrees with that found in the context of AdS/CFT correspondence. 

The paper is structured as follows: in Section 2 we present the general framework and analyse the 
field equations and boundary conditions exactly. In Section 3 we address the question of AAdS black hole 
solutions with complex scalar hair. In AD asymptotically flat space, the boson stars do not allow for black 
hole horizon inside them |35| . We generalize this result to include a negative A, for any D > 3. In Section 
4 the field equations are solved numerically. It is shown how the dimensionality of spacetime affects the 
properties of solutions and several cases are examined in detail. In section 5 the stability of the solutions 
is considered. The mass of the solutions and the boundary stress tensor are computed in Section 6. We 
conclude with Section 7 where the results are compiled. 



2 GENERAL FRAMEWORK AND EQUATIONS OF MOTION 

2.1 Basic ansatz 

Evolution of a minimally coupled complex scalar field $ in D dimensions is described by the action (through- 
out this paper we will use units in which c = Ti — 1) 



S D = -I d D X [^j^{R-2k) + L*)--^— I d D ^ x V~hK, (1) 

where the second term is the Hawking-Gibbons surface term I^Hj, K is the trace of the extrinsic curvature 
for the boundary dM. and h is the induced metric of the boundary. Of course, this term does not affect the 
equations of motion but is important when discussing quantum properties of gravitating solutions. 
In the lagrangian density of a complex self-gravitating scalar field reads 

£• - -V=gZ Or*-, +V(*)) , (2) 

where the asterisc denotes complex conjugate. In this paper we consider only the case V(&) = /i 2 $*$, 
where /i is the scalar field mass, without including a scalar self-interaction term. As found in |14| for four 
dimensional asymptotically flat solutions, although the inclusion of a A|$| 4 term drastically changes the 
value of the maximum mass and the corresponding critical central density of the boson star solutions, the 
qualitative features are essentially similar to the non-self interaction case. 

The Lagrangian density 10) is invariant under a global phase rotation $ — * <I>e~ M ; that implies the 
existence of a conserved current 

J* = igM , (3) 

and an associated conserved charge, namely, the number of scalar particles 

N = [ d^x^gEJ*. (4) 



The field equations are obtained by varying the action Q with respect to field variables gij and $ 

l>',j ; H ■ -V'/., = 8nG D Tij, (5) 

V 2 $-^ 2 $ = 0, (6) 
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where the energy momentum tensor is defined by 

Tij = 'l\'/t\, + **•*,< - '/,,i'/ i " ,< I , >'i , .». + /'"' •!' -)• (7) 

As usual with the scalar field, the Klein-Gordon equation is redundant. 

In this paper we assume a stationary ansatz $ = 0(r)e -1 " 4 describing a spherically symmetric bound 
state of the scalar field with frequency u>. Hereafter, a prime denotes the derivative with respect to r and 
an overdot d/dt. As we assume spherical symmetry it is convenient to use the metric form 

ds2 = W) + r2dnl - 2 F ^ e ~ 2S{r)dt2 > ^ 

where d£l 2 D _ 2 = cu a bdx a dx b is the line element on a unit (D — 2)-dimensional sphere and 

2m(r) 2Ar 2 , , 

F ^ = l -^- {D -2) { D-iy < 9 > 

The function m(r) is related to the local mass-energy density up to some D— dependent factor. Substituting 
in Q we find for the nonvanishing components of the energy-momentum tensor, 

for the energy density and the radial pressure while for the (D — 2) angular components we find 

T$ = P t = -F<f>» + ??f-<f?-f<l?. 

We remark that the radial component of the pressure P r and tangential components Pt are different, preclud- 
ing a macroscopic description of these configurations in terms of an effective perfect fluid energy-momentum 
tensor. 



2.2 Reduced action and virial relations 

If one is only interested in spherically symmetric solutions it is much simple to work with a reduced action 
where this symmetry of the spacetime is factored out. Expressing the curvature scalar R in terms of the 
metric function m(r) and S(r), we obtain the following expression for the effective action of our static 
spherically symmetric system 

S(m, 6, <j>) = I" ( { Jlz^l m 'e- 5 - r D - 2 e- s ( - ^£ + V ^)) ) dr , (1 ) 

J \ &TT& D \ F J J 

This form of the reduced action allow to derive an useful virial relation by using a scaling technique. In this 
way it is possible to better understand the reason for the existence of nontrivial solutions and to provide 
nonexistence theorems. 

We will use an approach proposed by Heusler in [37]. Let us assume the existence of a solution 
m{r),5{r),(j){r) with suitable boundary conditions at the origin and at infinity. Then each member of 
the 1-parameter family 

m\(r) =m(\r), 6\(r) = 8(\r), <j>\(r) = <f>(Xr) (11) 

assumes the same boundary values at r — and r — oo, and the action S\ = S[m\, 5\, (f>\] must have 
a critical point at A = 1, [dS/d\]\=i = 0. Thus we obtain the virial relation valid for a D— dimensional 
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spacetime 

* '"^ (<° " 3 » + ^ - (D-2KD-3) ) *" + < C - 1 >^> 

°-i- < " , ;f-. i '°; i ' tf ,)^=t (12) 

By using this relation we can also discuss whether the cosmological constant can support a real scalar field. 
In asymptotically flat spacetimes, a well-known result implies the absence of scalar solitons. It can be seen 
from l|12f) that in any dimension D > 3 there are no static, spherically symmetric scalar solutions with 
negative cosmological constant and a positive potential, independently of whether or not gravity is taken 
into account. 

Also, there are different proofs of the no-hair theorem for spherically symmetric scalar fields (see jSZj for 
a set of references). The above scaling argument does not exclude black hole solutions with complex scalar 
hair (to this end, one has to replace the lower boundary in the action integral l|l(Jf) by the horizon distance 
Th). However, by using another sort of argument, we will prove in Section 3 the absence of such solutions, 
which clearly indicates the limits of scaling techniques. 

For V{4>) = n 2 4> 2 , the relation (JUJ can be simplified by using the Klein-Gordon equation written in the 
form 

1/ D-2 -8T?fj,2\r\> D-2 ( tt±i2 , 2 ,2 J,2 a;2e25 \ /-, Q \ 

-[r e F(<p ) ) = r \F(p + fx <p - <p — — I ■ (13) 
Integrating the above relation and making use of the boundary conditions from Section (2.3), we find 

(f^ + ^ 2 -^^)=0, (14) 
which, together with 1|12[) gives the simple relation 

' * (d - + ^(p^^y - =*£^V) - » M 

In the above relation, the factor in front of <f) 2 has a negative sign. Thus, the existence of boson star solutions 
is supported by the factor (1 — m(3D — 7)/r D ~ 3 ) in front of cf>' 2 , which has no definite sign. 
The relation (|15|l particularized for D = 3 

/ 2 25a 2 \ 

dr re- 5 (1 - 2m)0' 2 + J ej? = 0. (16) 



dr r D - 2 -- 5 



F 2 

can be used to exclude the existence of three-dimensional asymptotically flat boson stars. 
2.3 Field equations and boundary conditions 

To perform numerical computations and order-of-magnitude estimations, it is useful to have a new set of 
dimensionless variables. Thus, we perform the rescalings r — ► r//i, <fi — * \J [^Gd ^' m ~~ * rn//i D ~ 3 and 
A — > A//i 2 , while the factor uj/fi is absorbed into the definition of the metric function S. 
Then we find the field equations 

D-2 / e 2S 6 2 \ 
m' = ^(V + 2 + ^j, (17) 

(e- 5 )' = r(e-^' 2 + e s ^), (18) 



{r D - 2 e- 5 F^)' = r D - 2 e- s d?[l- e —). ( 1 ! I ) 



2-> 
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Taking into account the explicit form for the metric and the scalar field we find the particle number 



47r (£>-i)/2 



N = Z,d-u I dr ^ rD ^ e -- ( 2 °) 



Similar to the A = case, we can define a star radius 

r / ,n_i T f , *7r I ,9 n_i e 



R = — I d D - 1 xrJ t J^B= — ^ / dr ^r - 1 — . (21) 

W V y T(^)N Jo V F 



Following the standard analysis, we can predict the boundary conditions and some general features of the 
finite-energy solutions. The boundary conditions for the system are the following. The non-singularity at the 
center of the star requires that 0(0) = ipo, 0'(O) = 0, m(0) = 0. We find also that the boundary conditions 
at infinity obtained for an asymptotically flat spacetime remain valid, with (f> = as the only acceptable 
value, and also that 

lim 0' 2 ~ 0(l/r D+2+e ), (22) 

r — >oo 

where e is a small positive quantity. 

The formal power series describing the above boundary conditions at r = is 

4>{r) = + ^ D "^ ] r 2 + 0(r% (23) 

m(r) = ^^ ) r D - 1 +0(r D ), (24) 

S(r) = 8o-\e 2S ^lr 2 + 0{t^). (25) 

The dimensionless value of the scalar field at the origin Lp$ is also used as a parametrization of the central 
density of the configuration. For A = 0, boson star solutions are characterized by an exponential decay of 
the scalar field, for which a mass term in the potential is responsible. The behaviour for A < is different: 
the analysis of the field equations as r — > oo gives 

4>{r) - 4> r c + . . . , 

m(r) „ M + (pc2 + Vfty c+D-i + , , , (26 ) 

V ' 2(2c+D-l) y ' 

where c = -| (d - 1 + ^(D - l) 2 +4/fc 2 ) , fc 2 = -2A/(D - 2)(D - 1) and ^ is a constant. Thus, the 
cosmological constant implies a complicated power decay at infinity, rather than an exponential one as found 
in an asymptotically flat space. 

We will prove in Section (6) that the constant M in the above relations 

is the ADM mass of a boson star up to a -D-dependent factor. However, in the discussion of numerical 
solutions we will refer to M to as the mass of the solutions. 



3 A NO-HAIR THEOREM 

In general, when a theory allows regular configurations, it also allows black hole solutions. However, this 
is not the case for the action discussed in this paper. In four dimensions with A = 0, Pena and Sudarsky 
have proven that there are no black hole analogues of the regular boson star configurations Therefore, 
the collapse of a boson star results in a trivial Schwarzschild black hole. A crucial point in their proof was 
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the asumption of asymptotic flatness. However, in the last years it has been realized that the introduction 
of a (negative) cosmological constant may allow for configurations forbidden in the asymptotically flat case 
(see, for example, the existence of dyonic regular and black hole solutions in EYM theory with A < 0). The 
dimensionality of spacetime may also affect the existence and the properties of solutions in a given theory. 
A nontrivial example is the gravitating SU(2) theory, which presents spherically symmetric particle-like 
solutions in D = 4 without counterparts in D — 3 j2Jj and D = 5 0Uj dimensions. 

In this section we prove the following straightforward generalization of the Pena-Sudarsky theorem for 
D > 3 and A < 0. 

Theorem: 

A D— dimensional, static, spherically symmetric, AAdS black hole spacetime, with regular event horizon, 
satisfying Einstein's equations with the matter fields fulfilling the weak energy condition (WEC) and its 
energy-momentum satisfying the condition 

Tl < T;, (27) 

is necesarily trivial (i.e. Tf = and the only black hole is Schwarzschild-AdS spacetime). 

Proof.: Since our arguments are very similar to |35j. we present only the main steps. From the (rr) 
and (tt) Einstein equations we can derive the following relations 

F' = i^T/-^ + ^(l-n (28) 
D~2 1 D-2 r K h y ' 

# = (W^F&-T r r ) (29) 
By using these relations and the conservation of the energy-momentum tensor T^ k — we find 



e s (e- s TZ)' = 



2Fr 



2Ar 2 

(D 3)(1 -F- — — — — )(T* - T;) + 2F(D - 2)(T 6 e ~ T r r ) 



(30) 



The WEC states that the energy density p = —T\ is semipositive definite and T^,Tg > T\. Thus, this 
condition and the assumption l|27l) implies that e~ s T^ is a nonincreasing function of r and e~ s T^(r) < 
e~ s Tr(rh). Here r% is the event horizon radius, corresponding to a zero of the function F. However, by 
expressing H30J) in terms of the proper radial distance dx = we find 

d 1 / 2Ar 2 \ 

e&-(eSTl) = -_ (D-3)(l-— — — — )(T*-T;) (31) 



dx y r) 2VFrY n {D-2)(D-3)' 



^ (2(D 2)(T e e T r r ) -{D- i)(Tl ~ T r r )) . (32) 

Since the left hand side of this equation is finite as r — > r%, we find T^(r^) = T t *(r/j) = — p(r^) < 0. However 
e^T^oo) = 0, so we conclude T, r (r) = 0. The WEC plus condition J2@ implies that T 6 B = Tf = also. 
The only black hole solution of Einstein field equations is thus the Schwarzschild-AdS spacetime. 

It can easily be seen that a complex scalar field with the lagrangian density J5J) and $ = §(r)e~ wt 
satisfies the conditions of the theorem above as long as V($) > 0. 

By using similar arguments we can also exclude the nonexistence of topological black holes. These are 
black holes with A < for which the topology of the horizon is an arbitrary genus Riemann surface, the 
(D — 2)-sphere being replaced by a (D — 2)-dimensional space of negative or vanishing curvature (sec 41 
for a review). 

Applying the same arguments in the case of a positive cosmological constant leaves open the question 
of the existence. We find only that the scalar field and the energy momentum tensor should vanish on the 
cosmological event horizon |42| while T£ is a strictly increasing quantity. 
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4 NUMERICAL SOLUTION 



While an analytical solutions to the coupled nonlinear equations i|17fl - l|19|) appears to be intractable for every 
dimensions, the resulting system has to be solved numerically. In this way we find that the boson stars with 
negative cosmological constant may exist in any dimension D > 3. 

A complete analysis of the complex correlation between the two parameters of the theory (ipo,A) for a 
given spacetime dimension is beyond the purposes of this paper. To compare numerically the results with 
those known in the asymptotically flat case we focused on solutions with D = 3, 4, 5 and have varied the 
parameter ipo for a set of A. 

For each choice of ipo, the eqs. (I17fl - I|19f) with the initial conditions ll'-'il) have a solution satisfying the 
boundary conditions (|26() only when e s (r = 0) takes on certain values. Different values of 5q corresponds 
to different number of nodes of the scalar field. We follow the usual approach and, by using a standard 
ordinary differential equation solver, we evaluate the initial conditions at r = 1CP 3 for global tolerance 
1CP 12 , adjusting for fixed shooting parameter and integrating towards r — > oo. We solve the equations and 
find the functions m, 8 and (j> and also the values of M, N, R. The energy of these solutions is always 
concentrated in a small region. 

Only the zero-node solutions will be given here. However, we found that similar to the D = 4, A = 
case <l> can have nodes, giving rise to excited states. 

4.1 Boson stars in four dimensions 

We start by discussing the better known case D = 4. The behaviour of M, Wasa function of (f>(0) is well 
known for the A = case; the mass and particle number rise with the increasing of </>(0) to a maximum 
of M max = 0.633 M Pl j \x and N max = 0.653 Af|, ; //i 2 . After this M and N drop, oscillate a bit and finally 
approach constant values independent of </>(0). 

From Figures lb, 4b it is clear that the introduction of A < does not change this behaviour qualitatively. 
The location of maximum shifts to higher values of </>(0) while the maximum mass and particle number 
decreases with the value of the cosmological constant. However, the general properties of the solutions are 
the same as for the A = case. The results for A = agrees with previous results in the literature. 

As a general feature, we have noticed a decresing of the maximal allowed value of the parameters N, M 
and a smaller ADM mass for the same central density as compared with A = case. For example, when 
A = —0.05 we have found that M max = 0.480 and N max — 0.435, while for A = —1, these values are 
M max = 0.239, N max — 0.121. Different values for the shooting parameter Sq are found. 

From Figure 2b it is evident that the characteristic boson-star masses decrease with increasing |A|. We can 
qualitatively understand this fact in the following way. Adding a negative cosmological constant corresponds 
to adding an attractive force. The Heisenberg uncertainty principles requires that a quantum state confined 
into a region of characteristic radius R has a typical boson momentum p ~ 1/R. Since for a marginally 
relativistic boson star p ~ /i we find R ~ l//x. The order of magnitude of the critical mass for the formation 
of a black hole can be estimated by comparing R with the Schwarzschild radius r s . Since the relation between 
Schwarzschild black hole mass and event horizon radius is 2M — r/i(l — A/SrfJ^ 1 ~ r/i(l + A/3r 2 ), we find 
a critical mass M c ~ (1 + A/3/i 2 )Mp l / fx, always smaller than the asymptotically flat value. 

Note that, again, for a finite value of the central density, the binding energy of the configurations Eb = 
M — N becomes positive, signaling the existence of configurations with excess energy. This property express 
a global instability against dispersion of particles to infinity. The excess energy is translated into kinetic 
energy of the free particles at infinity. 

4.2 Boson stars in five dimensions 

We integrated the equations of motion for a number of dimensions greater than four. We find that the 
general picture we presented for D — 4 remains always valid for higher dimensional spherically symmetric 
bosons stars. 
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The parameter M rises quickly with the initial value of the scalar field 0(0) at a maximum, afterwards 
drops, and has several local minima and maxima until it approaches an asymptotic value for large 0(0). A 
similar behaviour is found for the particle number. The location of the first peak again shifts to higher values 
of 0(0) for a negative A. Curiously, the five-dimensional asymptotically flat solutions present this peak for 
very small values of ifQ. 

The results for three diferent values of 0(0) and a varying A are presented in Figures 2c, 3c. As expected, 
the values of (M, N) for a given ip depends on the value of A. For a nonzero A it is necesary to establish 
new (and smaller) limiting values for the maximal mass and maximal particle number. For example, an 
asymptotically flat five-dimensional solution has M max = 3.803 and N max = 3.792; when A = —0.01 we 
have found that M max — 1.764 and N max = 1.667, while for A = —0.1, these values are M max = 0.893, 
N m ax = 0.719. We can heuristically understand this fact by using the arguments presented for D = 4. 
Typical results of the numerical integration are presented in Figures lc,4c. 

Note that again for a finite value of the central density the binding energy of the configurations E = 
M — Nfi becomes positive, signaling the existence of configurations with excess energy. 

4.3 Boson stars in three dimensions 

The (2+1) dimensional boson stars are a rather special case. Three dimensional gravity provided us with 
many important clues about higher dimensional physics. It help that this theory with a negative cosmological 
constant A has non-trivial solutions, such as the BTZ black-hole spacetime 02], which provide important 
testing ground for quantum gravity and AdS/CFT correspondence. Many other types of 3-D regular and 
black hole solutions with a negative cosmological constant have also been found by coupling matter fields 
to gravity in different way. (2+l)-dimensional stars with negative cosmological constant have been studied 
in ^102]. We notice also the existence of interesting BTZ-like solutions in a class of 3D gravity models in 
which the cosmological constant is induced |46j . 

Exact boson star solutions with A < (in the limit of large self interaction) are also known to exist j22 • 

The line element (JSJ reads in this case 

ds 2 = + r 2 d9 2 _ (1 _ 2m(r) _ Ar 2 }e -28 dt 2^ (33) 

1 — 2m(r) — At 

where the angular variable 9 is assumed to vary between and 2-7r. At the first sight, it may be that the 
boundary condition l|24|l is too strong and we have the freedom to add a constant to the function m(r) 
32 . However, if we assume that this metric posseses a symmetry axis (located at r = 0), and has no 
conical singularities, we have to impose the condition lim r _>o m(r) = 0. Asymptotically, the line element 
(|33|l approaches the BTZ metric 

rlr 2 

ds * = _J1 + r 2 d9 2 _ ( _ M _ Ar 2 )e- 25 dt 2 , (34) 

-A/ - Ar 2 

with M = — 1 + 2m(oo). We will find in Section 6 that the interpretation of the function m(r) and of its 
limit as r — * oo is more subtle . 

We have integrated the field equations with the initial conditions l|23|l for a number of A < and have 
found a different picture as compared to the case D > 3. The solution with ipo = corresponds to the global 
AdS^. For a given A, we find nontrivial solutions up to a maximal value of ipo, where the numerical iteration 
diverges. A divergent result is obtained also in the limit A — ► 0. The results of the numerical integration 
both for M and N and three values of A are presented in Figures la, 4a. In Figures 2a, 3a we present the 
variation of M and N with A for fixed values of the central density. 

We noticed that, for every A, 0o, the asymptotic value of the function m(r) is smaller that 1/2. However, 
the maximal allowed value of the particle number can be higher than 1/2 (for small enough values of |A|). 

Further details on these solutions as well as rotating 3D configurations will be presented elsewhere. 
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5 STABILITY OF SOLUTIONS VIA LINEAR PERTURBA- 
TION THEORY 



An important physical question when discussing selfgravitating configurations is whether these solutions are 
stable. For a boson star, the discussion of gravitational equilibrium is more complicated since we cannot 
apply the well-known stability theorems for fluid stars. When the spacetime is not asymptotically flat the 
stability analysis can also be a quite involved and subtle problem. 

A number of authors discussed the dynamical stability of four dimensional A = boson stars by means 
of linear perturbation theory [HJ-UHl- This lead to an eingenvalue problem, which is of Sturm-Liouville type 
and which determines the normal modes of the radial oscillations and their eingenvalues \ 2 . The sign of 
the lowest eingenvalue Xq is crucial; if Xo > then the star is stable, whereas it is unstable. In particular 
Gleiser and Watkins used a powerful method to discuse the stability of asymptotically flat boson stars, by 
studying the frequency spectrum of a class of radial perturbations j^2] . They have shown that the transition 
from stability to instability occurs always at the critical points of mass (or charge) against the value of the 
scalar field at the origin. A more powerful method based on catastrophe theory was also proposed 24 ; , with 
similar conclusions. 

Since the general picture we find for D > 4 is similar to the well known asymptotically flat solutions, we 
expect to find similar results for the stability question. Therefore we conjecture that M max represents the 
boundary between stable and unstable gravitational equilibrium. 

The method proposed by Gleiser and Watkins [221 will be used here to discuss the stability of D di- 
mensional solutions with A < 0. Since our approach is practically similar to that presented in |22j . we 
will present here the main steps only, emphasizing the points where the different asymptotic structure of 
spacetime is essential. We find that, for any dimension, the perturbation equations reduce to a set of two 
coupled differential equations. 

We consider the situation where the equilibrium configuration is perturbed in a way such that the 
spherical symmetry is still preserved. These perturbations will give rise to motions in the radial direction. 

To conform with previous works in D = 4, we use in this section the following parametrization for the 
spherically symmetric line element 



where v and A are functions of r and t only. 

The equations governing the small perturbations are obtained by expanding all functions to first order 
and then by liniarizing the equations We write therefore 



Similar to the A = four dimensional case, it is possible to reduce the full system of equations to only two 
coupled equations in two unknowns. If we choose to obtain these equations for 6<f>, it is convenient to introduce 
the notations 5<j){r,t) ~ f\{r,t) + icf>o(r)g(r,t) (in some cases, the analytically equivalent set of equations for 
<5A, fi may be more suited in the numerical computation). For the equilibrium configurations the metric 
functions are time independent, e A ° = 1/F, e v ° = e~ 2S F and the scalar field is <fi(r,t) — (f>o{r)e~ luJt . 

For this scalar field ansatz, the linearized (r, t) Einstein equation can be trivially integrated in time giving 



ds 2 = e x dr 2 + r 2 dfl 2 D _ 2 - e v dt 2 , 



(35) 



\{r,t) = \ (r) + 5\(r,t), 
v(r,t) = v (r) + Sv(r,t), 

where SX(r,t), 5v{r,t) are related to the perturbations of m, 5. The scalar field is written as 



(36) 




(37) 



6X = 



D - 2 



(38) 



The (r, r) and (t, t) Einstein equations give two more equations for the metric perturbations 




[(D - 3)e- Xo SX + re- x " 5\' - re- Xo \' SX] = 8ttG d e^duu 2 ^ 



(39) 
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+e- Xo S\cf > l - 2e~ v °uct> Q {uh - o g) - ^ X °^f[ - ^ 2 Mx 



ti\ 2r 
5v' = (D - 3)— + v' Q 8X + -p^f™G D e Xo 



- Vo u?<gbv + 2e- 1/ °w0 o (w/i - 0o5) 



><j>' 2 SX + 2e- x ° ( f>' fi-2f J L 2 cl> f 1 



(40) 



while the linearized Klein-Gordon equation for fx and g are 

f" + fx (^T^ + \H - A o)) - e X °-"°fx + he^i^e-^ - fi 2 ) - 2 we A °-^0 o5 
+^(8u' - SX') + e x °- va w 2 (j)v{5x - Sv) - ^V°0 O <5A = 0, 



g" + g' 



^—^ + l(v'a-K) + ^] -ge X °- ya +2ue Xo - vo ^- + \ue x °-' J( >{5 x -5 V ) = 0. 
r 2 0o J <Po 2 



(41) 



(42) 



The procedure now is simple. Using eq. 1)38(1 . it is possible to eliminate 5X from Q41J1 . An expression for 8v 
is obtained for eq. 14211 . We can also differentiate eq. I|42|l with respect to r and eliminate Sv' by using eq. 
(I40|l . In this way we find 



/" + f {^A + \H K)) he X °- Va + Ci(r)/i - 2^g" - C s (r)g' 

-ug' (^20 o + M^y^ + "'o- Ao)) = 0, 
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where 



o 



,\ -»o g l + 2e Xo-u ./ _ C2(r)g / _ _^ C3 ( r)/l = 0) 
00 0n 



(43) 



(44) 



G 2 



\ / 9 „ 9 . . . lQirGn ,,o,2D — 6 . , ,. 647rGn \ , 9 



= (D-2) 



D-2 
(K-"o) D-2 1 



D-2 



9 0o - A o) + 2(T i )- 2 T 1 
2 0o 0o 

16ttG d 



(45) 



+ ^(A - ^o)K - K + 4) - ^^r e A -^ 2 2 K - A + :), 



2L> - 6 S 



2 

w(20 o 



L> - 2 



>o' ^D-2 

167rGn , ,..2D — 6 . ,,. 327rGn ,, 9 % 
-r0 2 o ( + v' - X' ) - -—^^V" 



D-2 



We consider only radial perturbations, which conserve the total particle number N, as given by l|2(J|) . from 
which we can compute SN. Remarkably, we find that always the integrant of SN is a total derivative with 
respect to r 



The appropriate boundary conditions for r — * 00 are fx -> 0, g'r D+2c -> while /1 =const., r 2 5 ' 
origin. 

Now we assume that all perturbations are of the form e iX * , where x is the characteristic frequency to be 
determined. The self-adjoint system of coupled equation 143|) and (|44ll can be written in the compact form 



(46) 



at the 



with fi = g' and 



Lij — 



L i]f] = X 2 Mijfj, i,j = 1,2 



-i-Gxl + GxCx 2u±Gx4>, + GxG 
A. 

1 dr 



(47) 



dr 

-2LjGxfo£ + GxC 3 -£G 2 £ + G 2 C 2 



Mij = e Xo ~ Uo 



Gx 
G 2 
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where d = r D-2 e (v -\ )/2 ^ q 2 _ r £>-202 e 3( Wo -A o )/2 j which i s the required "pulsation equation". This 
sytem defines a characteristic value problem for x 2 which must be real (note the close analogy with the 
D = 4, A = case). With the above boundary conditions, it can be easily verified that both L y - and My 
satisfy the condition 



and similarly for My. The equations i|43[) and l|44|) can also be derived from a variational principle. Eigen- 
functions related to different eigenvalues satisfy also an orthogonality relation. The eigenvalues \ 2 form 
a discrete sequence Xo — Xi ^ X2 • • •• Thus, if the fundamental mode is stable (xo — 0)j then all radial 
modes are stable. Also, if the boson star is radially unstable, the fastest growing instability will be the 
fundamental one. We are interested in finding the values of </>o(0) for which the equilibrium configurations 
are stable. As extensive discussed in j22] . we can find some insight in this problem by considering static per- 
turbation. If for a given equilibrium configuration a radial mode changes its stability property (therefore its 
corresponding eingenvalues Xq passes through zero), then there will exist infinitesimally nearby equilibrium 
configurations with the same total mass and total particle number. Hence, if x 2 goes through zero we have 

dM/# (0) = dN/dMO) = 0. 

Also, if we begin with an equilibrium configuration with <pQ, the perturbed fields will describe another 
equilibrium configuration with 0o+<50o- However, the perturbed configurations must have the same charge as 
the equibrium configurations. Therefore, zero frequency perturbations will exist if and only if there exist two 
neighboring configurations with the same charge. This occurs at the extremum of N, where dN/d(j)o{G) = 0. 
This is a very general argument and aplies for every dimension. 

For all discussed cases except D — 3, the behaviour of the mass and particle number, parametrized 
by ipo suggest that the solutions with ipo smaller than ^ c (0) corresponding to the maximum mass, are 
stable, whereas the ones with bigger values than </> c (0) are unstable. Thus we expect that the value of ipo 
corresponding to the maximal mass is the boundary between stable and unstable equilibrium configurations. 

Of course, a complete proof still requires an analysis of the eingenvalues of the pulsation equations for 
given (D,A) in order to establish if they are real and positive for central densities smaller than C (O). 
However, it is a tedious procedure to find the characteristic values of x 2 f° r a given value of 0o (see 07] for 
a detailed discussion of this problem in D = 4, A = case). For our case, we solved the equations (Effl 
in a neighbourhood of C (O), for D = 4 and several negative values of A (we consider nodeless unperturbed 
solutions only). The boundary conditions are as follows: for r-»0we find 



We can use the linearity of the equations i|43[l . (|44|l to set /i(0) = 1. Dimensionless quantities are obtained 
by using the same rescaling discussed in Section 2. We then have two parameters left % 2 and q\ and the 
solution of the problem can be found by using a standard shooting procedure. The numerical results in this 
case confirm that x 2 > f° r central values of the scalar field smaller than </> c (0) while other configurations 
are unstable (% 2 < 0). In Figure 5 we present the result of the numerical integration for the perturbation 
of the metric function A for three diferent configurations with maximal mass (x 2 = 0). The variation of x 2 
with the central value of the scalar field is presented in Figure 6 for two negative values of the cosmological 
constant. We have all reasons to expect to find a similar result for any A < and D > 3. 



6 THE MASS AND THE BOUNDARY STRESS TENSOR 



An important problem of AdS space concerns the definition of mass and angular momentum of AAdS 
spacetimes. The generalization of Komar's formula in this case is not straightforward and requires the 
further subtraction of a background configuration in order to render a finite result. This problem was 
addressed for the first time in the eighties, with different approaches (see for instance Ref. |48l 149) ). 




h=Po+P2r 2 +0(r 3 ), g' = qi r + 0(r 2 ), 



(48) 
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At spatial infinity, the line element JSJ can be written as 



ds 2 = dsl + h^dx^, (49) 

where are deviations from the background AdS metric ds^. Similar to the asymptotically flat case, one 
expects the values of mass and other conserved quantities to be encoded in the functions h^ v . 

Using the Hamiltonian formalism, Henneaux and Teitelboim |3S] have computed the mass of an AAdS 
spacetime in the following way. They showed that the Hamiltonian must be supplemented by surface terms in 
order to be consistent with the equations of motion. These surface terms yield conserved charges associated 
with the Killing vectors of an AAdS geometry. The total energy is the charge associated with the Killing 
vector d/dt. A similar result has been obtained by Abott and Deser 0E|> m a more general context (see Ref. 
50 for a discussion of the relation between the Hamiltonian mass and the Abott-Deser one). This method 
requires to choose a reference background with suitable matching conditions. However, in some interesting 
cases (such as configurations with NUT charge) this choice is ambiguous. 

Another formalism to define conserved charges in AAdS spacetimes was proposed in [5Tj and uses con- 
formal techniques to construct conserved quantities. This construction makes no reference to an action. 
However, it yields the results obtained by Hamiltonian methods. 

A different approach has been proposed recently for locally AdS4 asymptotically geometries and gener- 
alized in for asymptotically AdS2« spacetimes (see also [EI])- It consists in adding an Euler (boundary- 
)term to the gravitational action |52| . The purpose is to guarantee that the action reaches an extremum 
for solutions which are locally AdS at the boundary. The resulting action gives the mass and angular mo- 
mentum as Noether charges associated to the asymptotic Killing vectors d/dt and d/dip. In this formalism 
the specification of a reference background is also unnecessary. The efficiency of this approach has been 
demonstrated in a broad range of examples, includding Kerr-AdS and Reissner-Nordstrom-AdS black holes, 
black strings and Taub-NUT and Taub-Bolt solutions. 

As expected, these different methods yield the same total mass for the spherically symmetric AAdS boson 
star solutions considered in this paper (for computations using Hamiltonian methods we consider the obvious 
vacuum AdSp background) 

M ADM = V ' - M, (50) 

07TO- D 

where n D „ 2 = 2^ D -^I 2 /T((D - l)/2) is the area of a unit (D — 2)-dimensional sphere. 

A procedure leading (for odd dimensions) to a different result has been proposed by Balasubramanian and 
Kraus [3^. This technique was inspired by AdS/CFT correspondence and consists also in adding suitable 
counterterms I c t to the action in order to ensure the finiteness of the stress tensor derived by the quasilocal 
energy definition |55| . These counterterms are built up with curvature invariants of a boundary dM. (which 
is sent to infinity after the integration) and thus obviously they do not alter the bulk equations of motion. 

Given the potential relevance of these boson star solutions in an AdS/CFT context, we present here a 
computation of the boundary stress tensor and of the total mass, by using the counterterm prescription. 

The following counterterms are sufficient to cancel divergences in a pure gravity theory for D < 6, with 
several exceptions (see e. g. |55]1: 



8ttG 



d D - x xy 

dM 



D-2 I 

1 1Z 

I + 2(D - 3) 



(51) 



Here 1Z is the Ricci scalar for the boundary metric h. Our conventions with respect to indices will be: 
{A, B, . . .} indicate the intrinsic coordinates of the boundary and {a, 6, . . .} indicate the angular coordinates. 
In this section we will define also I 2 = -(D - 2)(D - 1)/(2A). 

Using these counterterms one can construct a divergence-free stress tensor from the total action 
I=hu\k+Isw:t+Ict, by defining (see e.g. Ref.[57|) 

2 SI 1 D 2 I 

TAB = V^hSh^ - fagC*** - Kh ™ — kAB + —l EABh (52) 
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where Eab is the Einstein tensor of the intrinsic metric Kab ■ 

We can use this result to assign a mass to an AAdS geometry by writing the boundary metric in an ADM 
form 

h AB dx A dx B = -Nldt 2 + cr ab (dx a + N%dt)(dx b + N b dt) (53) 

and the definition 

M = I d°- 1 x^N^e. (54) 

Here e = u^u^T^ is the proper energy density while u M is a timelike unit normal to S. 

The metric restricted to the boundary \iab diverges due to an infinite conformal factor r 2 /l 2 . The 
background metric upon which the dual field theory resides is 

I 2 

j AB = lim -^h A s- (55) 

r — >oo T 

For the asymptotically AdSo solutions considered here, the (D — 1) dimensional boundary is the Einstein 
universe, with the line element 

jABdx A dx B = -dt 2 + l 2 dfl 2 D _ 2 . (56) 

If there are matter fields on A4, additional counterterms may be needed to regulate the action. However, 
we find that for a boson star spacetime in D — 3, 4, 5 dimensions, the prescription H51J) removes all divergences 
(this is a consequence of the asymptotic behavior of the scalar field (|26Jl ). The results we find for the boundary 
stress tensor of a boson star at large r by using the asymptotic expressions are 



T$ = l(M-±)5 ab + ..., (57) 

.(3) 



r;r = \{m-\)+..., (58) 



in three dimensions and 



Tif = ^% + -. (59) 
& = ^ + (60) 



for D = 4. The boundary stress tensor for the five-dimensional case is 

-,(5) -~ - I 2 - 1 

ab 

<5) 



T% = (M+-)-u ab + ... (61) 
3M 31 



By using the relation (|54|) . we can find the mass of boson star solutions, according to Balasubramanian and 
Kraus. The mass of a four-dimensional boson star is M — M, predicted also by (|50() . The result for D = 3 
is M = M — 1/2, while in five dimensions we find M = Madm + 37rZ 2 /32. However, in D — 5, the standard 
value is Madm- The additional constant 37W 2 /32 is the mass of pure global AdS$ and is usually interpreted 
as the energy dual to the Casimir energy of the CFT defined on a four dimensional Einstein universe |34| . 

For D = 3, we find that the mass of solutions computed in this way, is always larger than —1, with the 
extreme value M = — 1 corresponding to the global AdS^ space. Contrary to what happens for BTZ black 
holes, we do not find a mass gap between M = — 1 and M = 0. To compute the mass for this type of 
configuration, usually one uses a matching procedure on a surface separating the regions where the internal 
and external geometries are defined (see e.g. jS2])- The external geometry is taken to be the BTZ black 
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hole which gives a value for the mass of the solutions. The counterterm method gives, however, a rigurous 
prescription to compute the total mass of a 3D star, without the rather unnatural matching procedure. 

In light of the AdS/CFT correspondence, Balasubramanian and Kraus have interpreted Eq. I|52|) as 
< t ab >= 77=^ s-yAB 1 wnere < tAB > 1S ^ ne expectation value of the CFT stress tensor. Then, the 
divergences which appear are simply the standard ultraviolet divergences of a quantum held theory and we 
can cancel them by adding local counterterms to the action. Corresponding to the boundary metric lfST)|) . 
the stress-energy tensor tab for the dual theory can be calculated using the following relation 



However, further progress in this direction seems difficult, since we do not know the underlying boundary 
CFT for the action Q. In particular, we do not know what the gravitating complex scalar held corresponds 
to in CFT language. 



Many particle theories predict that weakly interacting bosonic particles are abundant in the universe, and 
that dark matter can be bosonic in nature. 

In this paper we investigated the basic properties of boson star solutions in the presence of a negative 
cosmological constant. Both analytical and numerical arguments have been presented for the existence of 
nontrivial solutions. A general virial relation has been found and particular cases have been discussed. 

The results we find for D > 4 are broadly similar to those valid for the known asymptotically flat, four 
dimensional case and there are only some small differences. Nontrivial solutions are again found for every 
value of the central scalar field and discrete values of frequencies while the scalar field vanishes at infinity. 
However, a nonzero A term in the action implies a complicated power decay of the fields at infinity, rather 
than exponentially as expected. Also the parameters of the solutions in the asymptotically flat case (mass, 
particle number, effective radius) are no longer valid and take new values for every choice of A. The mass 
of the star decreases with the cosmological constant strength |A|. Therefore the existence of a negative 
cosmological constant implies a decrease in the maximal mass of the star. In the sense of the no-scalar- 
hair theorem we have shown that, for this system, the only spherically symmetric black hole solution in 
D— dimensions is the Schwarzschild-AdS solution. 

These results are not so surprising, since the different asymptotic structure of spacetime obtained for a 
A ^ does not manifest here directly. The existence of the term V($)y/— g in the action (JIJ still implies 
$ — * asymptotically (since g rr gtt — > — 1), independently on the presence of A. 

A similar mechanism will act also for a Higgs scalar interacting with a nonabelian field in an AAdS 
geometry [SJ[8]. Here we consider only the better known case of a SU(2) field in four dimensions. In this 
case, a spherically symmetric purely magnetic YM field can be written as (see e.g. 0]) 



(where u> is the YM potential and are the usual Pauli matrices), which implies an YM effective lagrangian 
Lym — l/2e~ 5 (Fu' 2 + V(cj)), with V(uj) = (lj 2 — l) 2 /2r 2 . The scalar field effective lagrangian is L$ = 
r 2 e~ s (F(j)' 2 + V{(f))), with V(4>) the standard Higgs potential. The lagrangian of the system contains also an 
interacting term Li = e~ s T(uj)<j) 2 , where T(ui) = (ui + l) 2 for a doublet Higgs field |S] and T(u) = ui 2 for a 
Higgs field in the adjoint representation |2J. Again the finite energy requiement forces V{4>) — + as r — > oo, 
which implies that the Higgs field approaches asymptotically the vacuum expectation value. Moreover, the 
interacting term T{uj)4> 2 implies that the allowed values of the YM potential at infinity are similar to the 
A = case. In particular, the solutions are found again for a discrete set of the parameters that specify the 
initial conditions at the origin. Without a scalar field, the YM behavior is very different. The resulting YM 
equation reads {e^ 5 {Flo'))' — e~ s uj(oj 2 — l)/r 2 and the YM potential u> takes arbitrary values at infinity 
u) = loq + uti/r + 0(l/r 2 ). The solutions are found for continuous sets of shooting parameters and have a 
nonzero magnetic charge 01 . 




(63) 



r— ►oo 



7 CONCLUDING REMARKS 



A = - Lu(r)Tid9 + (cot6*T3 + uj(r)T2) sin 9 dtp 



(64) 
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Thus, it seems that when studying a gravitating scalar system, the solutions exhibit a generic behavior 
for A < 0. Given the presence of a cosmological event horizon, the case A > needs a separate analysis. 

We noticed, however, that there is a curious change of behaviour for D = 3, with somewhat different 
solution properties. 

We expect also to obtain a very similar qualitative behavior of the solutions when discussing a number 
of possible extensions of this theory, e.g. including a dilaton or a sclfmteraction term. 
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Figure Captions 

Figure 1: 

The mass-parameter M (dotted line) and the particle number N (solid line) are represented as a function of 
0(0) for typical boson stars with A < in three- (figure la), four (figure lb) and five dimensions (figure lc). 

In figures 1-4 the mass parameter M is given in units - — ^ — — ^ ^E'_ 3 > while the particle number N is 

given in units ^l^ (Ma)^, 
6 4r(u-i)/2) v ^ ; 

Figure 2: 

The mass-parameter M is represented as a function of the negative cosmological constant for different values 
of the central scalar field ip for boson star solutions in three- (figure 2a), four (figure 2b) and five dimensions 
(figure 2c). 



Figure 3: 

The particle number N is represented as a function of the negative cosmological constant for different values 
of the central scalar field tpo for boson star solutions in three- (figure 3a), four (figure 3b) and five dimensions 
(figure 3c). 

Figure 4: 

The particle number N is represented as a function of the effective radius R (in units 1/fi) for different 
values of the cosmological constant in three- (figure 4a), four (figure 4b) and five dimensions (figure 4c). 

Figure 5: 

The perturbation of the metric function SX is shown as a function of the radial coordinate for the funda- 
mental mode and several values of A. 



Figure 6: 

The squared frequency \ 2 of the perturbations is plotted as a function of scalar field at the origin for D = 4 
asymptotically AdS boson star solutions with A = —0.05, — 0.5. Note that % 2 is zero when $(0) ~ 0.48 for 
A = —0.05 and $(0) ~ 0.652 for A = —0.5, which corresponds to boson stars with maximum possible mass. 
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Figure lb. 
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Figure 2c. 
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